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SUMMARY 


FAFRAP is a computer program which calculates far-field radiated sound 
pressure levels from quantities computed by a NASTRAN direct frequency 
response analysis of an arbitrarily shaped structure. Fluid loading on the 
structure can be computed directly by NASTRAN or an added-mass approximation 
to fluid loading on the structure can be used. Output from FAFRAP includes 
tables of radiated sound pressure levels and several types of graphic output. 
FAFRAP results for monopole and dipole sources compare closely with an expli- 
cit calculation of the radiated sound pressure level for those sources. 


INTRODUCTION 


FAFRAP computes far-field radiated sound pressure levels using the 
Helmholtz exterior integral equation by numerically integrating fluid pres- 
sures and normal velocities over the fluid-structure interface of a finite 
element model. The numerical integration requires the XYZ coordinates, unit 
normal vector, tributary area, fluid pressure, and outward normal velocity for 
every grid point on the fluid-structure interface. ALTER statements in a NAS- 
TRAN direct frequency response analysis are used to obtain these quantities. 
Fluid pressures at the fluid-structure interface are computed directly by NAS- 
TRAN if an explicit fluid finite element mesh is used. Alternatively, FAFRAP 
will calculate fluid pressures if an added-mass approximation to fluid loading 
is used. 


THEORY 


Consider the arbitrarily shaped body in figure 1. Let z be the position 
vector to an exterior fluid point P, and z = ! zi . Let x be the position vec- 
tor to a point on the fluid-structure interface (with x = !xj), let r s z - x 
(with r = ij“! ) , an^ let £ be the unit outward normal at the location £. The 
time-harmonic (e ) pressure at z is given by the Helmholtz integral (ref. 1) 

p(z) = / g[iop v n (jc) + (ik + 1/r)p(j()cosS ](e - ^ r /4ir r ) dS (1) 
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Figure 1. Geometry for Far-Field Radiated Sound Pressure Level Computations 

where v n (_x) and p(x_) are the complex normal surface velocity and complex sur- 
face pressure, respectively, and k = u)/c, where w r circular frequency, and 
c = speed of sound in the fluid. 

Equation (1) can be simplified if only far-field locations are of 
interest. As \z\ + 00 , ik + 1/r ->• ik, and from the law of cosines, 
r z - x cosa . Therefore, at far-field locations 

p(z) = (ik e“ lkz /4Trz) /_ [pcv (x) + p(x) cosB ]e ikx COSa dS (2) 

o n — 

where cosB = (z/!z!)* n. 

One convention for presenting far-field pressures is as "sound pressure 
level (RMS) in dB relative to 1 MPa at 1 yard." Sound pressure levels (SPL) 
due to an excitation force applied as amplitude rather than RMS is obtained 
from eq. (2) by substituting z = 36 inches and by multiplying by the conver- 
sion factor 1 psi = 6.895x10 MPa to convert pressure p(^) from pounds per 
square inch (psi) to micropascals (uPa). Therefore, for \z\ = 36 inches 

SPL = 20 log ( ( 6 . 895 x 1 0 9 p(z))/^2) ( 3 ) 


FLUID LOADING 


Fluid pressures on the fluid-structure interface can be computed by NAS- 
TRAN if an explicit fluid finite element mesh is used. The finite element 
method that models the exterior surrounding fluid out to a predetermined dis- 
tance is described by Everstine (refs. 2,3). This method uses, as the funda- 
mental unknowns, the structural displacements and a velocity potential in the 
fluid. The outer boundary is terminated with nonreflective (wave-absorbing) 
boundary conditions, which assume that the outgoing waves are locally planar. 
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This approach to fluid loading results in an accurate model of the fluid mass 
at the expense of a much larger model due to the increased number of degrees 
of freedom introduced in modeling the fluid region. 

An alternative to using an explicit fluid finite element mesh is to use 
an appropriate added-mass approximation to fluid loading. The added-mass is 
applied to the grid points on the fluid-structure interface. For example, at 
low frequencies for a conical section, the effect of the fluid pressure is 
that of a mass-like impedance. Junger and Feit (ref. 4) show this impedance 
to be 


z = -iwn (4) 

a 

where the effective added mass per unit area is 

m a = (5) 


where p^. = fluid mass density, R = radius of conical section, and n = circum- 
ferential harmonic number. 


IMPLEMENTATION 


The numerical integration of eq. (2) requires, for each grid point on the 
fluid-structure interface, the XYZ coordinates, unit normal vector, tributary 
area, fluid pressure, and outward normal velocity. All of these quantities 
can be obtained directly from NASTRAN using the 0UTPUT2 utility module. The 
following ALTER statements will output the required data blocks on to the NAS- 
TRAN UT1 file. 

$ FAFRAP ALTER STATEMENTS, RIGID FORMAT 8, APR 84 VERSION 

ALTER 21 ,21 $ 

GP3 GE0M3 .EQEXIN .GE0M2/SLT .GPTT/S , N , NOGRAV/NEVER= 1 $ 

ALTER 55 $ 

SSG 1 SLT , BGP DT , CSTM , SIL , EST , MPT , GPTT , EDT , , C ASECC , DIT/ 

PG/LUSET/NSKIP 

SDR2 CASECC , CSTM , MPT , DIT , EQEXIN , SIL , GPTT , EDT , BGP DT , , , , 

EST , XYCDB , PG/0PG1 , , , , , /*STATICS*/S , N , N0S0RT2/-1 / 

S, N .STRNFLG $ 

ALTER 137 $ 

0UTPUT2 PG.BGPDT, EQEXIN ,FRL,UPVC $ • 

ENDALTER $ 

These ALTER statements allow a static unit pressure load to be applied to 
the structure during the dynamic frequency response analysis. The unit pres- 
sure load is applied to the fluid-structure interface of the finite element 
model. The components of the load vector created by the static pressure load 
are used by FAFRAP to compute the unit normal vector and tributary area of the 
grid points on the fluid-structure interface. The following is a list of the 
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quantities in the data blocks written with the 0UTPUT2 statement. 

PG - load vector components 
BGPDT - XYZ coordinates 

EQEXIN - internal to external grid point numbering equivalencing 
FRL - frequency response list 
UPVC - grid point displacements 

The displacements are converted to velocities by the relationship 
v : in u, where u is displacement and v is velocity. If an explicit fluid 
finite element mesh has been used, then the pressure at a fluid grid point on 
the fluid-structure interface is evaluated as the time derivative of the velo- 
city potential (ref. 2). If an added-mass approach to fluid loading is used, 
FAFRAP calculates pressure from the displacement. 

Several user-defined input parameters to FAFRAP control the number of 
far-field locations at which to calculate an SPL and the different types of 
output . 


OUTPUT 


Several types of output are available from FAFRAP. There are tables of 
computed values and three types of graphics output. Table 1 lists the SPL at 
far-field locations. The headings COLAT and LON refer to colatitudinal and 
longitudinal far-field locations, respectively. These tables are printed for 
each subcase and frequency. Table 2 lists phasor sum, RMS velocity, maximum 
SPL and where it occurs, maximum SPL in a horizontal plane and where it 
occurs, and radiated power for each subcase. Equations 6 and 7 define the 
phasor sum and RMS velocity, respectively, 


phasor sum = ?(v .A.) / A 

l n , l i 

RMS velocity = ( I v . ! ~k . ) / A 

l n , l l 


( 6 ) 

(7) 


where the summation is for all i grid points on the fluid-structure interface. 
Radiated power represents a summation of all pressure intensities in the far- 
field . 

A separate plotting program, FAFPLOT, was written to display SPL's in any 
of the three principal planes for any subcase or frequency. Figure 2 is an 
example of a polar plot of SPL generated by FAFPLOT. The two numbers in the 
lower left-hand corner refer to the subcase and frequency defined for that 
plot. The polar plot is useful in evaluating the sound pressure pattern gen- 
erated by the structure. 

Log-log plots can be generated for plots of pressure, velocity, or 
impedance at a grid point on the fluid-structure interface versus frequency. 

An example of this type of plot is shown in figure 3. The log-log plots are 
useful in evaluating the response of specific points on the fluid-structure 
interface for different load cases. 
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PATRAN (ref. 5) can also be used to display all of the SPL's in the far- 
field for one subcase and frequency as a color contour plot (fig. 4). This 
type of plot gives a good view of the overall radiated sound pressure pattern. 


COMPARISON TO ANALYTICAL SOLUTION 


Analytical solutions exist for the pressure fields produced by two simple 
radiators, the monopole and dipole sources. The equations defining the pres- 
sure fields generated by these sources were used to validate the results of 
FAFRAP. A NASTRAN analysis was performed for each of the sources to provide 
the necessary input for FAFRAP. The equations defining the pressure fields 
for a monopole and dipole source can be found in equations 4.15 and 4.75 
respectively, of Ross (ref. 6). After converting the equations to provide 
results in the correct units, the difference between the far-field radiated 
SPL calculated by FAFRAP and the values obtained from the equations was less 
than one percent. 


286 



REFERENCES 


1. Lamb H., Hydrodynamics , sixth edition, Dover Publications, New York, 1945, 
p. 498. 

2. Everstine, G.C., "A Symmetric Potential Formulation for Fluid-Structure 
Interaction," Journal of Sound and Vibration, Vol. 79, PP. 157-160, November 
1981. 

3. Everstine, G.C., F.M, Henderson, and R.R. Lipman, "Finite Element Predic- 
tion of Acoustic Scattering and Radiation From Submerged Elastic Structures," 
Twelfth NASTRAN User 's Colloquium , NASA CP-2328, National Aeronautics and 
Space Administration, Washington, D.C., pp. 192-209, May 1984. 

4. Junger, M.C. and D. Feit, Sound , Structures , and Their Interaction , The MIT 
Press, Cambridge, Massachussets , 1972, pp. 208,266. 

5. PATRAN User 's Guide , PDA Engineering Software Products Division, Santa 
Ana, California, 1984. 

6. Ross, D., Mechanics of Underwater Noise , Pergamon Press, New York, 1976, 

pp. 61 ,77. 


287 



6XAHPI 

? DF 

F A FRA P OUTPUT 


Ca/16/66 1 

2.29.57. 


soeciSF ;# 

FRt SU5NC Y - 

5 ' • ^ 

C jlAT 

LCN 

SPl ( C5 ) 

PHASE 

« >, 

• i 

92.651 

*i 2.3 

ifc.l 

• 0 

1:2* 22C 

89.1 

aC.C 

16.o 

102.622 

66.1 

16. C 

36.1 

1C3.39C 

92.5 

j . £ • » 

5 9.0 

1*. 2 • 9 •« 

97.7 

*e.» 

7<*t 

1 w 0 • t 2 7 

1C9.C 

ic.O 

9« • 

95.265 

112.1 

At • 0 

10 6.0 

99. C.7., 

126.9 

a 6 » c 

126.; 

91.223 

1 69 . 6 

ac. ( 

199.0 

91. C 27 

2C5.C 

*6 .1 

lit**. 

93*1 56 

2 31.3 

it* 1 

16-. .0 

89. £ £1 

2 53.7 

it,‘- 

196.0 

C4 , £ A3 

& * ♦ ** 

aC.C 

216. c 

92.516 

2 35.9 

iE*t. 

23*.o 

9C.697 

2 61.2 


252. 

82.266 

2 6-.-. 7 

iCl) 

Z7« .* 

69.366 

8b * 2 

a 6 . C 

26 6.’- 

92.926 

66.8 

At .0 

3o6.* 

97.929 

72.6 

aF.C 

321.0 

1 C C • 6 6 2 

76.6 

i E • V 

292. > 

1-2.267 

61 , 3 

j 6 • 0 

. 0 

1C5.197 

93.9 

3 6 . ■» 

*6 •*- 

1* 5.912 

99.. 

36«* 

26«« 

Iv 9.992 

55.3 

jc«C 

5 9 . v 

1 c 3 • e 8 9 

'•2.8 

36.0 

/ c . C 

::z.C72 

7c • 8 

3t,< 

9- • 

9<?.7t7 

c o • 2 

3t *C 

106.C 

96.222 

126.2 

J t • 0 

126.0 

86.6*9 

159.8 

3 t »' 

199.* 

99.7- 3 

17^.3 

3 fc.C 

162.0 

1 C t .222 

191. A 

36. C 

19 C. C 

1CC.C22 

196.2 

3 6.0 

198.0 

99.111 

i£9*. 

36. C 

216. C 

97.392 

1 90 . 7 

36.t. 

229.0 

99,2 92 

191.2 

36.* 

252.. 

6 6.9. 6 

1 53 .1 

36*0 

27C.0 

62.629 

55.8 

j6.o 

2c 6 •• 

89, £3' 

32.7 

36.0 

206. v 

99,726 

23.2 

j fc • t 

229.0 

lo2.582 

35.9 

36. C 

292.1 

1C 9. 3C6 

39.5 


Table 1. Sound Pressure Levels in the Far-Field 


288 



i' X AM PL t nf FAF9AP OUTPUT 


1 X 

/U/Cfc 

12.29.57. 










sc 

FRfcQ 

prIASC’R SP 


PMS VtL 

MAX spl 

cm at 

L 0*! 

H" 7 IP L 

COL M 

Lf r 

f r w 1 ,< 

i 

50. C 

l.l? 0 E-< a 

9 

, C 8 St - c. 

95.995 

199. 

1 9 e . 

91. 91.. 

3C 4 

9* . 

2 .:-!»• - 3 

X 

51 

1.2f26-< t 

5 

.2 3c F-Ot 

57.6S 6 

199. 

1 * 6 . 

5 2 • t 9 

3t 4 

5 6 . 

9,696 i — < 3 

X 

5 2.6 

1 .991 F— 1 6 

5 

.92i £- t 

C9.99C 

199, 

19 e. 

99.273 

364 

9 f . 

t ,111 \ -. 3 

i 

53.0 

1.770S-C5 

6 

. 7 0 1 ~ “ 0 1 

7 

2 9 9, 

i -i-3. 

96.1 ; 

Jo« 

R . 

9.5s 7-— > 

i 

5 9.6 

2.1 s^r. 

7 

• Sv.9t- 6 

I-, 3.2 '• 1 

199, 

i >e. 

97. ?c: 

!: o * 

90. 

1 . 9 5 •• P -i .* 

1 

5 5 • 

2 . 7 31 E-Ofc 

9 

• 9 25c--;.‘i 

lui .989 

19 9, 

1 9 8 * 

; 9 • K 7 

- c * 

9* , 

'•Hit- 2 

X 

5 6.. 

i.63ve-< i 

1 

.1 5Sf-f«r 

1C 7 . > 3 9 

199. 

1 9 6 . 

i r 2. o 7 

3 c • 

r n . 

j « 5 3 1 — 2 

X 

57. ( 

5.231 E-< C 

1 

.9916-. 5 

Hi .6 .9 

199, 

190. 

19.625 


9 ( . 

7 * 3 0 f' r — •. i 

i 

5F.0 

C.390F-C8 

2 

•0 93t— u5 

1 19.299 

.99. 

1 )8. 

■7.6.9 

i * 4 ♦ 

i 7< • 

.. » 5 ** . ; " ‘ i 

1 

59. v 

*•96 6t-< 7 

3 

» 0 1 3 E — : 

UP. f . 7 

199. 

198. 

1 1 i . Sit 

3 ^ • 

5 C • 

3.515 i- 1 

1 

fcf! • v 

1.62U-C7 

4 

• 266 E— 05 

121.721 

19 9. 

1 9 6 . 

i 1 1 , • 9 • 

**5 6 * 

5* . 

7 . 1 '■ 7 c -< 1 

i 

fci *o 

e.572f-('H 

3 

• f 51 v—06 

12C.121 

199. 

17 8. 

119. 953 

3 6* 

5v . 

3.95 3s— 1 

1 

62.1 

Z.P19E-I b 

? 

,t99c- 

R6 .£91 

199. 

196. 

il .953 

3 6 « 

SC. 

2.55 1*6— \ 1 

i 

53. C 

1. 3 77f-M: 

1 

. fc72t- r 5 

.'_9.958 

199, 

; 9 F . 

; 7 . 6 f • 


5« . 

1.9i 5:— i 

X 

6:4.1 

J .3 6 71—0 6 

i 

.9 2 6 1- 

12 2.659 

3 99. 

: 9p. 

• 6 . ; • ' 

.U . 

Si . 

5.75 < i . 

1 

65.0 

1.10*t-0b 

1 

•lilt-. 5 

111.29? 

299. 

19 6. 

i- 9. 5; 7 

1 1 • 

5 • . 

7 . i 5 1 — 2 

X 

66. 

i ,0 99 1 — C f- 

9 

,699 r -C6. 

11C. 210 

199. 

19 6. 

R i . 9 ‘ 'i 

3 1 • 

S 5 ft 

5 . r. 1 9 ! — l 

i 

6 7 . i. 

9. 6i.3L- : ’ 9 

8 

• BctE-t t 

32 i 

199. 

19 6. 

1- 2.529 

3t ♦ 

Sf c 4 

9 . 6. 2 6 : -v 2 

i 

fc f .0 

9.1634— 99 

7 

•3t iL-:t 

. t . 6 * ). 

i 9 9 . 

19 6. 

* . - 4 it 7 

1 r. • 

• 

c 

X 

69. t 

t .<! 7t 9 

5 

. 5 c J r — ■ 6 

)’ 7.969 

199. 

? 9 6. 

i ■ i • 2 - - 

X C 9 

Si 4 

2.5 * 6 1 -> ; 

X 

70.6 

C.CPV: -I'C 

* 

. 697^-06 

li 7.9 ;9 

199. 

1 7£. 

) • 4 H i 

c • 

<< 4 

j » ( 6 : — ■ ? 


Table 2. Summary of FAFRAP Results 
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Figure 3. Log-log Plot of Pressure Versus Frequency at a Grid Point 
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Figure 4. Contour Plot of Sound Pressure Levels in the Far-Field 



